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ABSTRACT 



Most, if not all, stars are now believed to produce energetic outflows during 
their formation. Yet, almost 20 years after the discovery of bipolar outflows from 
young stars, the origins of this violent phenomenon are not well understood. 
One of the difficulties of probing the outflow process, particularly in the case of 
massive embedded stars, is a deflcit of high spatial resolution observations. Here, 
we present sub-arcsecond-resolution mid-infrared images of one massive young 
stellar object, AFGL 2591, and its immediate surroundings. Our images, at 
11.7, 12.5 and 18.0 //.m, reveal a knot of emission Ri6" SW of the star, which may 
be evidence for a recent ejection event or an embedded companion star. This 
knot is roughly coincident with a previously seen near-infrared reflection nebula 
and a radio source, and hes within the known large-scale CO outflow. We also 
flnd a new faint NW source which may be another embedded lower-luminosity 
star. The IRAS mid-infrared spectrum of AFGL 2591 shows a large silicate 
absorption feature at 10//m, implying that the primary source is surrounded 
by an optically thick dusty envelope. We discuss the interrelationship of these 
phenomena and suggest that mid-infrared imaging and spectroscopy provide 
powerful tools for probing massive star birth. 



Subject headings: stars:pre-main-sequence-ISM: jets and outflows- 
stars: individual: AFGL 2591-stars: circumstellar matter 
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1. Introduction 

It is now generally agreed that most, if not all, stars produce energetic outflows during 
their formation (see Bachiller 1996 for a review). Since the discovery of outflows almost 20 
years ago (Snell, Loren & Plambeck 1980), nearly 200 examples have been found, most of 
them associated with low-luminosity young stellar objects (YSOs). Outflows can be studied 
using a variety of tracers including high- velocity CO line emission (Bally & Lada 1983), 
shock-excited H2 line emission (e.g.. Lane 1989), and optically visible Herbig-Haro (HH) 
jets (Reipurth 1989, 1999). 

All indications are that flows emerge bipolarly from a stellar or circumstcllar region. 
However, the physical processes giving rise to the outflows are not well understood. One 
suggestion is that a primary wind originating at the star (e.g., Shu et al. 1991) or at the 
accretion disk (e.g., Pudritz & Norman 1983) is responsible. In both scenarios, the fast 
well-coUimated wind sweeps up ambient molecular gas in its vicinity, forming two cavities in 
opposite directions from the star. The displaced molecular gas constitutes the CO outflow 
and ionized "knots" of gas manifest themselves as HH objects. Outflows probably occur 
while the star is still accreting material, and may play a significant role in ejecting material 
of excess angular momentum. 

The presence of "bullets" in the CO emission and shocked "knots" visible in the optical 
and near-infrared suggests that mass ejection from YSOs often occurs in bursts rather 
than at a steady rate. Such ejection events appear to be intermittent, and the interval 
between two successive events is on the order of 10^ years. It is intriguing that the masses 
and timescales of bullets are similar to those of outbursts observed in FU Ori stars. The 
FU Ori eruptions are thought to be due to large increases in the accretion rate through a 
circumstcllar disk (Hartmann & Kenyon 1985). If the two phenomena are related, it may 
be possible to probe the accretion history of protostars using detailed observations of their 



-4- 



outflows (e.g., Reipurth 1989). 

One of the barriers to better understanding the outflow process is the hmited spatial 
resolution of millimeter-wavelength CO line observations. The situation is particularly 
dire in the case of massive young stellar objects in which large extinction {A^ ^ 10 mag) 
prevents high-resolution optical imaging. 

Here we present sub-arcsecond-resolution mid-infrared images of one high-luminosity 
(Class O/I) YSO, AFGL 2591, and its immediate surroundings. AFGL 2591 was first 
recognized as an outflow source by Bally & Lada (1983). The CO outflow is larger than 
1' in extent, with the blue lobe extending west or southwest and the red lobe in the 
northeast direction (Lada et al. 1984; Mitchell, Hasegawa, & Schella 1992; Hasegawa & 
Mitchell 1995). Optical and near- infrared observations have revealed molecular hydrogen 
emission (40" east and west, Tamura & Yamashita 1992) and HH objects (10" and 20" 
west /southwest, Poetzel et al. 1992) near AFGL 2591. The distance to AFGL 2591 is not 
well established. Assuming a distance of 1 kpc, it has a luminosity of ~ 2 x and an 

estimated stellar mass of 10 Mq. 



2. Observations and data reduction 

AFGL 2591 was observed on June 3 and 4, 1999 using MIRAC3, the newest version of 
the University of Arizona/SAO mid-infrared camera (Hoffmann et al. 1998) at the 3.0-meter 
NASA Infrared Telescope Facility (IRTF). MIRAC3 uses a Boeing HF16 128x128 Si: As 
blocked impurity band detector. On the IRTF, MIRAC3 has a plate scale of 0.33" /pixel, 
providing a total held of view of 42" x42". This pixel scale ensures Nyquist samphng of the 
diffraction-limited point spread function (PSF). 

We obtained images of AFGL 2591 in MIRAC3 11.7, 12.5 and 18.0 //m 10% passband 
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filters, with a total on-source integration time of 600, 370 and 540 seconds respectively. The 
standard stars a Her and /3 Peg, observed before and after the science target, were used for 
flux and PSF calibration. 

We used a standard nodding and chopping technique to remove the background signal, 
dithering the source on the array to obtain sub- pixel sampling of the PSF. The chop 
frequency was set to 5 Hz, with a throw of 20" in the N-S direction. The nod throw was 
also set to 20" but in the E-W direction, in order to have all four chop-nod beams inside 
the field of view of the array. Each individual nod cycle required a 10-second on-source 
integration, and the procedure was repeated for as many cycles as needed to obtain the 
requested total integration time. 

To reduce the data, we have developed our own software written in C and IDL. We 
first subtracted the chop-on from the chop-off frames for both nodding beams. The two 
images thus obtained were then substracted one with the other, in order to get a single 
frame in which the source appears in all four beams (two negative and two positive). We 
then applied a gain matrix, derived from images of the dome (high intensity, uniform 
background) and the sky (low intensity, uniform background), to fiat field the chop-nodded 
image. 

This procedure was repeated for each of the nodding cycles for which the source 
was observed. A final high S/N cumulative image was then obtained coadding together 
all beams, each recentered and shifted on the source centroid. This last coadding was 
performed on a sub-pixel grid having the size of one-fifth of the original MIRAC3 pixels, 
thus providing a final pixel scale of 0.066" /pixel. A mask file to block out the effects of bad 
pixels and field vignetting was also created and applied in this stage, preventing individual 
rejected pixels from contributing to the final image. 

The same observing and reduction procedure was also used for the reference stars, to 
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ensure a uniform treatment of the source and the standards. 

3. Results and Discussion 

Figure 1 shows the final MIRAC3 images of AFGL 2591. The central source is not 
clearly resolved. It has FWHM of 0.84", 0.87", and 1.21" at 11.7, 12.5 and 18.0 /xm 
respectively. The FWHM of the reference point source j3 Peg in the same filters is 0.82", 
0.85" and 1.11" respectively. An extended knot of emission Ri6" SW from the central source 
is visible in all three images. A compact faint source is seen Rill" NW of AFGL 2591. 
Figure 2 is a schematic drawing of AFGL 2591, not to scale, with the various observed 
components indicated. 

In Table 1, we list mid-infrared fluxes for the three sources as well as the position 
angle of SW and NW sources (counter-clockwise from North) and their separation from the 
central source. The photometry was measured inside apertures of (diameter) 18" (total flux 
of main and SW sources), 9" (SW source alone) and 5" (NW source). The zero magnitude 
flux densities for the three filters are 29.2, 25.7 and 12.6 Jy respectively. We estimate an 
error of 10% for the photometry, mainly due to the uncertainty of the reference fluxes and 
the difficulty of isolating the sources when doing aperture photometry. Lower fluxes of the 
main source were found by Lada et al. (1984) by using broad-band filters at 11.4 and 12.6 
fj,m. The discrepancy (40% and 8% in the two filters) can be explained by the presence 
of a broad absorption feature in the source spectrum at 10 /xm, which affects broad-band 
photometry more than our 10% filters. By fitting the 12.5 and 18.0 /xm fiuxcs with a black 
body, we derive color temperatures of 300-350K for the central source, and 100-120K for 
the SW knot and the NW source. 

Persi et al. (1995) presented the only previously published mid-infrared image of 
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AFGL 2591. Their observations at 11.7 /im were carried out on the 2.1-meter S. Pedro 
Martir telescope in Baja Cahfornia, Mexico with an image scale of 1.34" /pixel, and did not 
reveal the SW knot or the NW compact source seen in our MIRAC3 images, due to the 
smaller field of view of their detector. 

3.1. Central source 

The fiux of the central source rises dramatically from 11.7 to 12.5 //m, implying that its 
11.7 iim fiux may be suppressed by a large silicate absorption feature at ~10/im. Indeed, 
the low-resolution mid-infrared spectrum obtained by the InfraRed Astronomy Satellite 
(IRAS) shows deep absorption at that wavelength. Thus, it is hkely that the protostar is 
surrounded by an optically thick dusty envelope, perhaps in free-fall collapse onto the star 
(van der Tak et al. 1999). Our observations limit the radius of the envelope to ;^3" (or ^ 
3000 AU, assuming a distance of 1 kpc). 

Figure 3 shows the IRAS spectrum of the source, and the best fit wc obtain by 
modeling the circumstellar dust envelope with the radiative transfer code DUSTY (Ivezic et 
al, 1999). The fit shown is for a spherically symmetric envelope with a radial density profile 
p (X around a central black body source at effective temperature T^ff = 25,000 K (van 
der Tek et al. 1999). We adopted the silicate dust opacity by Ossenkopf et al. (1992), which 
provides a good fit to the 10 iim sihcate feature. The fit for p oc r~^'^ (i.e., constant infall) 
is not as good, implying that the infalling envelope is likely to be surrounded by a fair 
amount of non-infalling cloud material. The model spectrum is not particularly sensitive to 
the Tgff of the central source. Our best fit to the IRAS spectrum suggests an Ay on the 
order of 100 (derived from the best fit optical depth). 
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3.2. NW source 

AFGL 2591 is generally considered to be an example of relatively isolated massive star 
formation. However, our detection of a faint point-like source 11" NW of it suggests there 
may be other embedded low-mass stars in its vicinity. A preliminary look at Two-Micron 
All-Sky Survey (2MASS) J,H,K images reveals at least two more point-like sources within 
40" of AFGL 2591. Deep mid-infrared images of a larger area are needed to search for a 
possible cluster of low-mass protostars in this region. 

3.3. SW knot 

The SW knot we detect appears close to the inner portion of a loop first observed 
in the near- infrared by Forrest & Shure (1986). The near- infrared emission is assumed to 
arise from scattering of photons from AFGL 2591 off dust grains at the surface of a cavity 
cleared by the outfiow. The loop, also seen in NH3 observations by Torrcllcs ct al. (1989), 
is indeed located in the blue outflow lobe, i.e., the one directed toward us. Using Draine & 
Lee (1984) opacity at 18/xm and assuming T =100K and d —1 kpc, we estimate the mass of 
this knot to be >O.2M0. 

The SW knot is also coincident with a radio source as well as a knot of H2 emission. 
The radio source, detected using interferometers at the Very Large Array (Campbell 1984) 
and the Owens Valley Radio Observatory (van der Tak et al. 1999), has an approximately 
flat spectrum between 6.1 and 0.3 cm, and is interpreted as free- free emission from an 
optically thin HIT region. The H2 emission is seen in an image taken with the BEAR 
imaging spectrometer on the Canada-Francc-Hawaii telescope using a narrow band filter on 
the 1-0 S(l) line, after the continuum is subtracted (G. Mitchell, personal communication). 
The existence of the H2 line implies shock coUisional heating of the molecular gas in 
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the outflow, which provide the necessary excitation energy {T ^ 2000 K, see e.g. Smith 
& Brand, 1990). On the other hand, the dust responsible for mid-infrared emission is 
characterized by a much lower local equilibrium temperature of ~ 100 K, as inferred 
from our photometry. It is likely that the H2 emission comes from the shock itself, while 
the mid-infrared color temperature is an average through the dust clump. 

Previous near-infrared observations have also revealed a second, larger loop of emission 
at a projected separation of ~19" from the star (Burns et al. 1989; Minchin et al. 1991; 
Tamura ct al. 1991; Poetzel, Mundt, & Ray 1992). The presence of multiple shells/loops is 
thought to be evidence for episodic mass outflows from AFGL 2591, as are molecular bullets 
seen in CO emission (Bachiller et al. 1991) and shocked knots observed in the optical and 
the near-infrared (Reipurth 1989) in other well-known outflow sources. Up to ~ IA^q 
could be ejected in each eruptive event, though a few tenths of a solar mass is more typical 
(Bachiller 1996). 

If we assume that the SW knot is associated with a recent mass ejection from AFGL 
2591, we can estimate its dynamical timescale. Following Minchin et al. (1991), we take 
that the outflow is inchned from the plane of the sky by 55 degrees, and estimate ~ 10^ 
AU as the true distance from the star to the SW knot. Given a true outflow velocity of 26 
kms~^ (based on the line-of-sight velocity of CO gas observed by Bally & Lada 1983), then 
we get a dynamical timescale of ~2000 years. For the same assumptions, the outer loop 
seen in the near-infrared would have a dynamical timescale of ~6000 years. Such timescales 
are consistent with the estimated 10**-10^-year total duration of the energetic outflow phase 
of YSOs (Bally & Lada 1983). 

Another possibility is that the SW source itself is an embedded protostar surrounded 
by an ultracompact HII region. Some Herbig Ae/Be stars -LkHa 198 and LkHa 234, for 
example- are known to have heavily embedded mid-infrared companions associated with 
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an outflow, a radio source, and HH objects (Lagage et al. 1993; Cabrit et al. 1997). The 
SW source is well-resolved in our mid- infrared images, with a "coma" -like shape. So, if it is 
another protostar, its envelope may be affected by the outflow from AFGL 2591. 

In summary, our observations suggest that mid-infrared imaging and spectroscopy can 
be powerful tools for probing the immediate circumstellar environment of embedded stars 
driving outflows. Such outflows appear to be rather common around young stars and are 
very hkely associated with accretion processes from the early stages of protostars. They 
may play an important role in determining the mass of the star-disk system, and affect the 
evolution of the surrounding molecular cloud. Yet our understanding of the outflow process 
is far from satisfactory. High-resolution imaging with current mid-infrared cameras and 
future millimeter-wave interferometers could help develop a more unified picture of infall 
and outflow in star formation. 

We are grateful to George Mitchell, Charles Lada and Paolo Persi for helpful discussions 
and to an anonymous referee for sensible suggestions. Wc thank the IRTF staff for their 
outstanding support. The research was supported by the Smithsonian Institution and 
NASA grants to the SIRTF-IRAC project. 
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TABLE 1 
MID-INFRARED PHOTOMETRY 

Source 11.7/im (Jy) 12.5;um (Jy) 18.0/im (Jy) PA Separation 

A 440 (±44) 745 (±75) 753 (±75) — — 

B 18.4 (±2) 15.9 (±2) 79.6 (±8) 235° 6" 

C 1.8 (±0.2) 2.1 (±0.2) 16.2 (±2) 325° 11" 
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Figure Captions 



Figure 1. MIRAC3 images of AFGL 2591 at (a) 11.7//m, (b) 12.5yum, and (c) 18.0//m. 
Contours are plotted for 10 and 15 a levels in panels (a) and (b), and for 5, 10 and 15 a in 
panel (c). 

Figure 2. A schematic drawing of AFGL 2591 with the various observed components 
indicated. (Not to scale.) 

Figure 3. IRAS mid-infrared spectrum of AFGL 2591 {solid line) and the best-fit 
model {dashed line) for the 10/im silicate feature. The heavy absorption is presumably due 
to an optically thick dusty envelope around the massive central star. 
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